The sudden approximation in photoemission and beyond 
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The conditions behind the sudden approximation are critically examined. The fuzzy band expres- 
sion is derived in detail from first principles. We go beyond the sudden approximation to account 
for both extrinsic losses and interference effects. In an extension of earlier work we discuss both core 
and valence satellites including both intrinsic and extrinsic amplitudes, and high energy excitations 
as well as low energy electron-hole pairs. We show how the extrinsic losses in photoemission can 
be connected with the electron energy loss function. This is achieved by three approximations, to 
connect the dynamically screened potential in the bulk solid to the loss function, to account for 
the surface, and lastly to extrapolate zero momentum loss data to include dispersion. The extrinsic 
losses are found to be weak for small loss energies. For larger loss energies the extrinsic losses can 
be strong and have strong interference with the intrinsic losses depending on the nature of the solid. 
The transition from the adiabatic to the sudden regime is discussed for atoms and localized systems 
and compared with the situation for solids. We argue that in solids for photoelectron energies above 
some 10 eV the external losses are mainly connected with spacially extended excitations. We discuss 
strongly correlated quasi-two-dimensional solids with Bi2212 as example, and find an asymmetric 
broadening of the main peak from shake up of acoustic three-dimensional plasmons. In the super- 
conducting state the loss function is assumed to have a gap, which then leads to a peak-dip-hump 
structure in the photoemission spectrum. This structure is compared with a corresponding structure 
in tunneling. 



I. INTRODUCTION 

Photoemission spectroscopy (PES) has become an im- 
portant tool to analyze electronic structure. The accu- 
racy of the experimental energy resolution is now in the 
meV range, and sample control has reached a high level 
of reliability. Yet the theory used to interpret these high 
quality experimental results has great weaknesses. Often 
only the spectral function is compared to experiment. 
This means that e.g. the momentum dependence of the 
dipole matrix elements is neglected, as well as the ex- 
trinsic losses. Use of the spectral function is fine if we 
have well-defined sharp quasi-particle peaks, and want to 
do band structure mapping. As the experimental reso- 
lution increases and line shapes are resolved this level of 
approximation becomes increasingly doubtful. That is in 
particular true for strongly correlated systems like high 
temperature superconductors (HTC) which have impor- 
tant features that cannot be described by quasi-particle 
peaks. 

Presently PES for HTC materials is often compared 
to spectral functions for models like the t-J one for a 
two-dimensional (2D) copper oxide layer. Such a model 
can only (at best) describe a narrow low energy region 
of the total excitation spectrum. The omitted degrees of 
freedom describe higher energies, and also collective low 
energy modes with three-dimensional character. As far 
as we know, there has been no systematic attempt to de- 
velop a theory that both gives a good description of the 



strongly correlated system and also takes account of the 
remaining correlations. This is an important and urgent 
problem if we want to have a true a priori understanding 
of e.g. HTC materials. In lack of a systematic theory 
we sketch below an approach to this problem built on 
physical intuition. 

The 2D problem is normally treated with a model 
Hamiltonian based on a small number of orbitals {4>2d}- 
The parameters in the model Hamiltonian can be ob- 
tained from " down- folding" the basis functions in a band- 
structure calculation by using e.g. developments of the 
very efficient LMTO formalism (TJ. The interesting cor- 
relation effects come from the degeneracy in the prob- 
lem, i.e. the large number of possibilities to occupy the 
different mi and m s states of the localized muffin tin or- 
bitals connected with the down-folded part of the band- 
structure. There are also efficient methods based on the 
LMTO formalism to handle a wide energy range which 
treats valence and conduction bands on an equal footing 
so that the wave functions are orthogonal [^J. Let's call 
the additional set of states beyond {</>2d} for {4>3d}- We 
assume that {</>3.d} is strictly orthogonal to {</>2d} despite 
any complications that can come from down- folding. The 
one electron states {4>2d} give rise to a huge number of 
configurations which then all are orthogonal to the con- 
figurations that can be built from the set {4>sd}- We 
also assume that we can write the state vectors as a di- 
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both the 2D and 3D parts are correlated state vectors. A 
parametric dependence of the 3D state on the 2D state 
is obvious when we have a core hole in the 2D system, 
but is actually important also in general. 

When we develop our theory we will not make any spe- 
cific assumptions about the correlations in the 2D system. 
We will assume that the excitations in the 3D system are 
extended in space and can be described by a quasi-boson 
model. This is possible if we only want to describe PES 
for energies larger than about 10 eV where a strongly 
correlated localized system already has reached the sud- 
den limit, m The treatment here of strongly correlated 
systems is based on our recent work. W 



II. THE "ONE-STEP" THEORY 

Bandstructure and surface effects in photoemission are 
today often treated in the so called "one-step model" pi- 
oneered by Pendr y M , for recent reviews see e. g. 
Refs. 0, @, §, In the one-step model the im- 

portant dipole matrix elements are included, as well as 
surface photoemission (emission due to the gradient in 
the electromagnetic field, for recent work see e. g. Ref. 
PI), but not the extrinsic losses. 

A complete theory was developed by Caroli et al in 
1973 |l2| based on the three current expression of Schaich 
and Ashcroft. [|l3| The Caroli et al theory uses the 
Keldysh formalism |lj] for non-equilibrium Green's func- 
tions. It is exact and one can e g identify the diagrams 
which give extrinsic losses. Caroli et al did not however 
extract any practical calculation scheme from their the- 
ory. 

Important steps towards a practical photoemission the- 
ory were taken by Feibelman and Eastman jl5| and by 
Pendry J5|, ||, who discussed the lowest non-trivial di- 
agram in the Caroli et al expansion, assuming the fi- 
nal states to be (time-inverted) LEED scattering states. 
Pendry and collaborators developed an explicit multiple 
scattering formulation and computer programs rooted in 
earlier work on the LEED problem. |r| A careful dis- 
cussion of the Caroli et al theory and the basics of pho- 
toemission was given by Almbladh 1985 leading to some 
explicit approximations. fl7| 

With todays precise experimental data on photoemis- 
sion the need for a comprehensive and practical theory is 
great. The one-step theory has serious limitations both in 
its theoretical foundation and in its scope. Thus Pendry, 
Feibelman and Eastman just assumed that the photoclec- 
tron states should be damped. That this is correct was 
shown later by Almbladh, however he uses the Keldysh 
diagram expansion which is rather involved. More se- 
rious is that the one-step theory neglects the extrinsic 
losses. These losses were early recognized to be impor- 
tant, and in core electron photoemission they are rou- 
tinely included. That they also are important in photoe- 
mission from an extended state (valence electrons) is not 



as widely recognized. 

A full theory for photoemission has a number of in- 
triguing features. We will display some of them by 
starting from the Golden Rule expression rather than 
the Keldysh expansion, and postulating a simple model 
Hamiltonian. Exact results can then be obtained by el- 
ementary means. The model Hamiltonian will be moti- 
vated a posteriori by discussing cases where the results 
can be obtained by other means. 

The Golden Rule expression for the photo-electron cur- 
rent can be obtained from scattering theory |Q (we use 
atomic units, e = h = m = 1, and thus e.g. energies are 
in Hartrees, 27.2 eV), 

J k (o;) = ^|(iV-l, S ;k|A|Ar)| 2 5( w - ek + £ 5 ). (1) 

s 

The final states \N — 1, s; k) have the proper boundary 
conditions for scattering states. All states are correlated. 
Further e-^ is the energy of the photo-electron, £k = k 2 /2 
p9[ , e s gives the energy of the final state of the solid as 
e s = E (N, 0) — E (N — 1, s), ui is the photon energy, and 
A the optical transition operator, 



\N) is the initial state. The exact expression for Jk (ui) 
is rather intractable, and we have to dev elop approxi- 
mations that can be handled. HJ, ||, HI For 

non- 
interacting electrons where \N — 1, s;k) and \N) are 
Slater determinants, it is easy to see that \N — 1, s; k) = 
c\x s \N), and thus Eq. [I] reduces to the well-known form 

occ 

J k ((j) = ^ l A k s | 2 S (UJ - £k + Es) , 

s 

where the index "s" now stands for an occupied one- 
electron state (e.g. a core electron state or a state in the 
valence band). 

The final state can be written as, 



\N-l,s;k) 
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(H — E) 



cl\N-l,s), 
(2) 



where H is the fully interacting Hamiltonian includ- 
ing target electrons and photoelectron, and E = lu + 
E (N, 0) = £ k + E (N - 1, s) . Further \N - 1, s) is the 
target state after the photo-electron has left, and cre- 
ates the photoelectron. The states corresponding to c£. 
are so far undefined, except that they are time-inverted 
LEED states with asymptotically a plane wave part 
exp(ik • r). Usually in scattering problems (H — E) is 
replaced by the interaction V between the scattering par- 
ticle and the target. This is not possible when the scat- 
tering particle is identical with particles in the target. 
One way to avoid this difficulty is to regard the scatter- 
ing particle as different, and do an anti-symmetrization 
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of the solution at the end of the analysis, which is the 
route followed in the standard treatise of scattering by 
Goldberger and Watson, (pi 

At high energies we can regard the scattering particle 
as different from the particles in the target since the ex- 
change coupling becomes negligible. We can then replace 
(H — E) by a potential V which describes the coupling 
between the photoelectron and the density fluctuations 
in the solid. Q This is true for all systems, also strongly 
correlated ones. For strongly correlated localized systems 
the interaction with the photoelectron becomes small al- 
ready at energies of about 10 eV. || As discussed in the 
introduction we will take the total wave function as a 



3D 
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For 



direct product, 1*404,31,52} = |*2D,si) 
energies above 10 eV we only need to consider the in- 
teraction between the photoelectron and the excitations 
"s 2 ". 

If the (H — E) term, which gives the extrinsic scatter- 
ing, is neglected, and if we assume that c k \ N) = 0, the 
amplitude becomes 

(N -l,s;k I A| N) = ]T (k I A| j) (N —l,s \ Cj \ N) . (3) 
j 

The photo current in Eq. [TJ is then given by the spectral 
function Aij (ek — to), 

J k H = J2 (k I A| i) Ay ( Ek - lj) (j |A| k) , (4) 



where 



N-l,s 



s\ Cj \N) 



xS (lj - s s ) . 



(5) 



Eq. |] is the well-known intrinsic approximation for the 
photocurrent. The photoelectron states |k) are time- 
inverted LEED states (with no magnetic field this just 
means complex conjugation). A LEED state has an in- 
coming and a set of scattered plane waves outside the 
solid, matched to a solution inside the solid with a po- 
tential Vh + S (ek), where Vh (r) is the Hartree poten- 
tial for a finite solid with a surface and E the complex 
selfenergy. In the next section we will show how the per- 
turbation V gives rise to extrinsic energy losses and the 
selfenergy term S (ek). 

In the one-electron approximation the spectral func- 
tion is diagonal and has just one quasi-particle peak, 
(lj) = S (lj — Ei) Sij. In most applications the spec- 
tral function is taken as diagonal. This is e g generally 
justified in the GW approximation. [^2| A core electron 
spectral function for a metal A c (lj) has an asymmetric 
peak (MND singularity ||^]) followed by a set of plasmon 
peaks. The asymmetry is due to shake-up of electron- 
hole pairs and phonons. [^4| For a conduction electron 
in a metal the spectral function A^ (lj) looks quite simi- 
lar with a broadened and asymmetric quasi-particle peak 



and a pronounced satellite structure. The broadening 
and asymmetry however go to zero at the Fermi level. 
Semiconductors have no low energy electron-hole excita- 
tions and thus lack the MND singularity, but still have 
pronounced plasmon satellites. Transition metals have 
about the same total strength in the satellite structure, 
but lack the pronounced high energy peaks of sp-metals 
and valence semiconductors. Strongly correlated materi- 
als may lack the quasi-particle peak altogether, and have 
all the spectral strength in the satellite. 

For core electron photoemission the current in 
the intrinsic approximation is simply Jk(w) = 
|(k |A| c)| 2 A c (ek — lj). For valence electrons we need a 
more involved discussion. We denote the initial state "j" 
in Eq. |^ with a vector 1 which also includes a bandindex, 
and Eq. [I] becomes 



JkH=^|<k|A|l)| 2 A 1 (e k -c). 



(6) 



For energies close to the quasi-particle energy E\, the 
spectral function A\ has the form (neglecting the asym- 
metry, which can be pronounced away from the Fermi 
level) 



Ai (lj) 



K (LJ- El) + r 



2 ' 



(7) 



which together with Eq. |^ gives the one-step approxima- 
tion. The photoelectron wavefunction |k) inside the solid 
is determined by the potential Vh (r) + E (r, r'; ek), where 
Vh is the Hartree or Coulomb potential. The real part of 
the non-local optical potential or selfenergy S is usually 
approximated by a local potential and combined with Vh 
to give an effective potential V e f / . Replacing — Im E by 

its expectation value Tf > 0, and setting |k) = k\ inside 

the solid, we have 



(t+V eff -iT f ) 



ek 



(8) 



where t is the kinetic energy, — — (k |ImE (ek)| ky > 

and ek = k 2 /2. The photoelectron wavefunction k 

inside the solid can be written as a combination of Bloch 
functions with complex k vectors, such that it correctly 
matches the plane waves outside the solid. If one Bloch 
wave dominates [ p5[ its complex momentum k follows 
from 



E. 



ek + iTf, 



where EL is the complex quasiparticle energy having both 

k 

a shift in its real part from Re E and an imaginary part 

T j. We write k=k|| + z (jt z + ikf \ , where k|| is parallel 

to the surface and z is a unit vector normal tojhe sur- 
face. Due to momentum conservation we have kn = k|| . 
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In the limit of small Tf we can expand £L 

k 

i Imk 9£L/<9k, and thus we have 



E, 



Rok 



dE„/dk z 



Mahan p6j and in more detail Spanjaard et al [£5| stud- 
ied the matrix element ^k|A|l^) in the limit when the 
damping is small. The integration parallel to the surface 
gives momentum conservation, k|| = 1||. The integration 
in the z direction gives a surface contribution plus contri- 
butions from the unit cells where there is bulk symmetry. 
With bulk symmetry all cells give the same contribution 
but for a factor, and we can sum to have 



klAll 



klAll 



(9) 



klAll 



mrf 



1 — exp 



(lz ~ 



k)c 



where c is the cell length in the z direction and ^k | A| 1^ 

is the integral over the first (closest to the surface) bulk 
unit cell. Band indices and 2D reciprocal lattice vectors 
have to be supplemented where appropriate. With small 

kj and neglecting the surface contribution /k |A| 1 

sharply peaked as function of l z around 

l° z = mod(fc z ,27r/c), 
and one can approximate 



klAll 



k|A|l° 



(l°-l z ) 2 +(kf) 



2 k,,,!,. 



where is a wavefunction in the occupied valence 
band, and 

1° = k,| + 21° = 1° (k) . 

The surface contribution to the photocurrent can be com- 
parable to the bulk one, and interference between the 
two terms in Eq. ^| is then important. jO], |^7j Collect- 
ing our results the contribution to the photocurrent from 
the quasi-particle peak is 



J k (uj) 



'A 



k I Al 1 



1 



(e k -Ei-lo 

where I\ = T\. 



With Ti small we can expand e k — E\ — uj = 
— (l z — ll) dE\jdl Zl where l\ is given by the resonance 
condition 



E, 



0. 



(10) 



We define 



l 1 =k|, +zll =1 : (k,c 



When also Ti is small we can integrate over l z keeping 
everything constant except the resonance denominators. 
We then have a convolution of two Lorentzians, which 
gives a new Lorentzian with the sum of the widths 



Jk (uj) 



(11) 




dEr 



dl z 



Z\ i Ti 



klAll 



where 



fcf = 



dE x /dl z 



i=i° 



k 1 - 



dE k /dk z 



kj 



k=Rck 



The momentum parallel to the surface is conserved, 
k|| = k|| , and thus immediately follows from experiment 
while the determination of k z requires knowledge of V e f / 
and thus e.g. ReE. There are however more involved 
experimental techniques to obtain also k z . |2£| As before 
band indices and 2D reciprocal lattice vectors have to be 
supplemented where appropriate. 

The transmission problem related to the matching of 
the photoelectron wavefunction inside the solid and the 
plane waves outside has been studied by Spanjaard et al. 
p5[ who conclude that for Ag the transmission factor is 
about 0.5 at all angles and for all energies of interest. 
The connection of the "fuzzy band" expression and dif- 
ferent modes of experiment has been discussed in detail 
by Smith et al. (g9[ The deviation from free electron 
energies for the higher bands, which shows up as the dif- 
ference between k and k, persist quite high up in energy. 



III. EXTRINSIC LOSSES AND THE SUDDEN 
APPROXIMATION 

The importance of extrinsic losses was recognized 
early, and the landmark paper here is the discussion by 
Berglund and Spicer in 1964 of the "three step model". 
pl| The three steps are photoexcitation, transport to 
the surface, and passage through the surface. In the 
last two steps the electrons can loose energy (extrinsic 
losses). Core electron photoemission data are often an- 
alyzed taking extrinsic losses into account. Methods to 
do this have been developed in a long series of papers by 
Tougaard and collaborators pS| . Strong extrinsic losses 
also occur in valence electron spectra. In Fig. 1 we show 
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the XPS valence band spectrum from Na metal found 
by Hochst, Steiner and Hufher |3^]. More details can be 
found in Refs. |3^] , fl35| |. It is remarkable that the struc- 
ture from the plasmon satellites is as strong as the main 
band (first peak). Actually the valence electron spectrum 
is similar to the core electron spectrum. This similarity 
is very close, the full drawn line in Fig. 1 is obtained by 
convoluting the Na core spectrum with a valence band 
spectrum obtained by weighing the partial densities of 
states (shown in the inset) with coefficients used as free 
parameters (s:p:d=l:5:9). This picture is confirmed by 
theoretical work. 1391, |B7f 




Binding energy LeV ) 



FIG. 1. XPS valence band spectrum of Na metal. Points 
are the experimental data and the full curve is the generated 
spectrum. The broken curve is the fitted background. The in- 
set shows the total and partial densities of states (DOS) which 
were used to obtain the generated spectrum. Full curve, total 
DOS; dotted curve, s partial DOS; broken curve, p partial 
DOS; chain curve, d partial DOS. [33] 

The plasmon satellites in the core spectrum are easy to 
understand, the sudden creation of a core hole changes 
the potential on the valence electrons and we have a 
"shake-up" effect. The surprising fact is that the sudden 
removal of a valence electron has a very similar shake-up 
effect, it is as if the valence electron hole were localized. 
A quantitative analysis of core electron satellites showed 
that both extrinsic and intrinsic losses contributed com- 
parable amounts. |3£| Since the valence satellites could be 
reproduced using the core spectra, we expect that also for 
valence electrons both intrinsic and extrinsic losses will 
contribute. 

A theoretical analysis of satellites in simple metals was 
done by Penn both for core spectra, |3!| and valence spec- 
tra, |56| . He solved a transport equation with the spectral 
function as source term. The spectral function for core 
electrons had been obtained by Langreth. Gof For valence 



electrons Penn obtained a spectral function quite similar 
to the core one by doing lowest order perturbation theory 
using an effective Hamiltonian based on the Bohm-Pincs 
plasmon theory. Penn's result clearly demonstrated that 
both intrinsic and extrinsic effects contributed to the 
satellites. 
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FIG. 2. Intensity of the first bulk plasmon satellite relative 
to the main peak for the Al 2p core electron spectrum. [45] 

As analyzed by Chang and Langreth plj , jf|| there 
must be interference between the intrinsic and extrin- 
sic processes since they lead to the same final state 
for the system. In particular they discussed "sp-solids" 
and found the interference effects to be strong and fade 
off only at kinetic energies in the keV region. A de- 
tailed and quantitative theory of interference effects in 
a core electron spectrum was developed by Inglesfield us- 
ing electron-plasmon coupling functions. [43], Q] Exper- 
imental results for the onset of the first p asmon satellite 
as function of kinetic energy have been obtained by Flod- 
strom et al. Q (see Fig. 2), and of others (see e. g. Ref. 
fS5f). Their results are well reproduced by Inglesficld's 
theory Q 

A more general approach based on a quasi-boson 
model, fluctuation potentials and the GW approxima- 
tion was developed by Hedin and collaborators. [M , 
|57j, J22j These results also show that both extrinsic 
and intrinsic effects are important, and that interfer- 
ence persists in simple metals to quite high kinetic en- 
ergies. This approach has a more rigorous basis and 
should be quantitatively more reliable. In particular it 
handles not only plasmon losses but also the low energy 
electron-hole losses which lead to the Mahan-Nozieres- 
deDominicis (MND) singular threshold behavior. |^3| We 
will now discuss photoemission using this quasi-boson 
model. 

Following earlier work by Hedin et al 20 1 , [^l| , |22| we 
use a quasi-boson model Hamiltonian, 



H = H o + V, H = H . 



syst 



(12) 



where 
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H. 



syst 



qW 



For the photoemission amplitude we now have (sup- 
pressing the strongly correlated part in the problem), 



E 



e k c k c k> v 



T. vq K 



v q 



kk' 



Here indices "I" label states below and "k" states above 
the Fermi level. Thus "I" can stand for a core electron 
state or an occupied Bloch state in the valence band, and 
"k" is a state above the Fermi level (but possibly below 
the vacuum level). States below the vacuum level (oc- 
cupied and unoccupied) are modified by the presence of 
the surface and have exponentially decaying amplitudes 
outside the solid. States above the vacuum level have 
the photoelectron boundary condition, i e they are time- 
inverted electron scattering states (LEED states) with 
free electron energies, £k = k 2 /2 and are plane waves at 
large distances from the solid. As basis set for the one- 
electron wavefunctions we take the eigenfunctions of the 
Hartree-Fock Hamiltonian. ]23] We will suppress spin 
variables. The term V describes the interaction between 
the photoelectron and charge density fluctuations (see 
below) in the solid, i. e. it is the interaction which causes 
the extrinsic losses and also the damping of the photo- 
electrons. Further a q stands for quasi-boson type excita- 
tions like electron-hole pairs, plasmons etc., and oj q is the 
quasi-boson energy. The quasi-boson parameters are re- 
lated to the dynamically screened potential W(r, r'; u>) as 
described in Appendix A. The potentials V q (r) and the 
wavefunctions can be chosen to be real functions, which 
we use in the derivations. The final expressions are valid 
also with complex functions. 

The model Hamiltonian we have defined clearly can 
only give an approximate description of a real system. 
Thus we obviously do some double counting of the 
Coulomb interactions. Despite this the predictions from 
the model Hamiltonian are basically sound. For self- 
energies it correctly describes the GW approximation, 
and it gives realistic results for both intrinsic and extrin- 
sic effects in core electron photoemission. p3] The intrin- 
sic satellites in valence electron photoemission for simple 
metals come out qualitatively correctly. fl37| Further the 
high energy limit in photoemission is correctly described. 
pl| Clearly a quasi-boson model cannot handle strongly 
correlated systems, and we will discuss these separately. 

Returning to Eq. || we take | N — 1, s) as an eigenfunc- 
tion of H syst with eigenvalues u s , and thus 4 \ N — 1, s) 
is an eigenfunction of Hq with eigenvalue Sk + lu s . Since 
E = £k + <jJ s we can replace (H — E) with V. The states 
| A — l,s) have a hole after the emitted photoelectron 
plus a number of boson type excitations. For core elec- 
tron photoemission the bosons have to be calculated in a 
potential with a core hole present, while the initial state 
| A) is the vacuum state of bosons with no core hole. Also 
in valence electron photoemission, as we have discussed 
earlier, the hole after the valence electron has a strong ef- 
fect on the quasi-bosons leading to a comparable strength 
in the satellite. 



(A-l,s;k|A|A) 



(13) 
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N-1,8 



Ck 



1 + V- 



E-H 
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The physical significance of this expression is quite clear. 
The operator c\cj lifts one electron from the state "j" 
into the state "i ". The "virtual" electron "i" propa- 
gates by (E — if) -1 and is then scattered by V into the 
state "k " which is measured by the detector. The state 
"i" must be a photoelectron state, which we denote by 
k'). Since we never have more than one photoelectron 
we can work with a product representation, 

4 | A - 1, s) = \N - 1, s) |k> , cl | A) = | A) |k'> . 



We can now rewrite Eq. Il3 as 



(A - 1, s; k |A| A) 
Wk (N-1,8 

Jk' x x 
x(k'|A|j) 

( k ( A — 1 . >• 



(14) 



1 + V 



1 + V 



1 



E- H + ir] 



N 



k' 



1 



E- H + irj 



N ) A 



noting that |k') (k'| acts as a deltafunction. 

The potential V in Eq. [l4| gives the two effects of a 
loss potential, the energy loss scattering and the related 
reduction in the intensity of the no-loss peak. The lat- 
ter comes through a partial summation over the no-loss 
intermediate state, which gives a damped photoemission 
wavefunction. Formally this summation can be done by 
use of the Feshbach projection operator technique. The 



result is (see e g Ref. 21 



(A-l, S ;k|A|A) 



(15) 



N —l.s 



1 + V 



1 



E - QHQ + it] 



A ) A 



k ) is the wavefunction discussed in connection 



where 

with Eq. ISl and Q is the projection operator Q = 1 — P, 



with P = \N — 1, s) {N — 1, s\. In principle 
self-energy, 



k) and the 



E = PVQ 



1 



E - QHQ + in 



-QVP, 



depend on the excitation state of the target, "s". To 
lowest order we have 14711 



S(r,r'; W )=^^ (r) 



1 



-h-< 



ii] 



-V (r') , (16) 



G 



provided the quasi-bosons in the state " s" are extended 
and thus a finite number of them have a negligible influ- 
ence. Also when we go to higher orders the result remains 
independent on the index " s" . The imaginary part of £ 
in Eq. [Hj] agrees with the imaginary part in the GW ap- 
proximation for the self-energy jl?]], while the real parts 
differ slightly. 

To estimate the extrinsic and intrinsic losses we ex- 
pand Eq. 15 to lowest non-trivial order in V q (see Ap- 



pendix B). To lowest order in V there is one contribu- 
tion from V (E — Hq)~ (extrinsic loss) and one from 
(N — 1, s \cj \ N) (intrinsic loss described by the spectral 
function). In both cases "s" is a state with one quasi- 
boson (q) and one hole. The result is given in Eqs. 
and |^|. The selfenergy £ in the denominator does not 
appear in the simplest lowest order perturbation theory 
but is obtained in a systematic expansion. jf7| 

We first discuss valence electron photoemission when 
the index I is a vector index 1. From the discussion at the 
end of Sect II we see that the matrix element | A| 1^) is 

peaked around 1 = 1°. The no- loss contribution to Jk (lo) 



j£(w) =£|(k|A|l) S(to-e k + e l 



klAir 



(l»-l z f + (k\) 



T S (lo - e k + £i) dl z 



k |A| 1° 



\de\ldl 7 \ 



(iv-iir + ikj 



which is the QP particle result in Eq. [Tl] in the limit 
when Ti (and Si) are infinitesimal. 

For the loss contribution to Jk (uj) we first consider 
only the intrinsic part (c.f. Eqs [l] and |4g| ), 



4 M = J2 



E 



-<k|A|l')ttf 



i'i 



LOq - £l + £\> 



8 (uj - io q - e k + £i) ■ 



(17) 



The quadratic expression in the fluctuation potentials 
can be written in terms of the imaginary part of the self- 
energy 

Im E M « M = vr V^V^S (u + u q - ei) , 



provided lo < fx, (see Eq. 19 in Ref. |22|). We can ex- 
press Im £ (lo) in terms of Im W (uj) as is clear from Eq. 
KS.This gives 



1 {k|A|l')(k|A|l")lm£,,,„( ek -^) 



IT 



(lo - e k + ey) (uj - e k + sv 



The condition uj < fi for Im E {oj) translates to e^—uj < fi 
in the above equation. This means that uj > e^—fi, where 
lo is the photon energy, which clearly is satisfied. Since 
Im Ei'i» {lo) is diagonal to a good approximation (Ref. 
||) we take 



J2» = £E- 

7T z — ' 



klAir 



Im Ei' (ek — lo) 



{uj - e k + El' ) 
We now make the following Ansatz for Jk (lo) 



(18) 



J k (w) = 



E 

i 



klAll 



^ t(U)-E k +gl)t 

2tt 



(19) 



^ /• Im Si (ei — lo') ( e 



1 W 



x exp 



The integrand in the exponent is the function Si (— t) 
discussed in Rcfs [El and , 



Si(t) 



which for large t is 



1 /■ Im£i(ei - w) (e lwt - lWw 



Si (i) = —T\ \t\ + iot\sgn (t) — n\. 

The quasiparticle behavior is set by the large t limit. Ne- 
glecting the asymmetry index ay we have 



i 



klAll 



r, 



(« - e k + EiY + r 



2 ' 



which agrees with Eqs |^ and M s ince the factor exp (— n\) 
equals Zi to leading order. ]37j] The Ansatz to leading 
order also reproduces the intrinsic loss J k (lo) in Eq. [l8| 
apart from a normalization factor. The Ansatz further 
agrees with the result for intrinsic losses obtained ear- 
lier. |57| To have the intrinsic core electron spectrum we 
just replace the index 1 by c in Eq. [l^, and we have no 
summation. 

For core electrons we can obtain both the intrinsic and 
extrinsic contributions [14|. a detailed derivation is given 
in Appendix A in Ref. \W. We found that Eq. 43 can be 
written 



(k) 



V q (r c ) 



fc Z) K|A[c 



1 f z ° 

+ -/ dzifc (z)V(q z ,Q,z)ij<(z) 
™ J -00 k * 

x («,K-Q|A|c), 

provided we make a number of simplifying approxima- 
tions. We have e. g. taken the crystal as translation- 
ally invariant parallel to the surface. K and Q are the 



7 



photoelectron and fluctuation potential momenta paral- 
lel to the surface, and r c is the position of the atom from 
which the core electron is excited. Further V (q z ,Q, z) is 
the fluctuation potential V q , and ip~ (z) and ?p< (z) are 

damped Blochfunctions, ip£ decreasing towards the sur- 
face, and ip> decreasing towards the inner of the crystal. 
The momenta k z and k are equal when the quasi-boson 
energy and momenta are zero. We also take the dipole 
matrix elements as equal, and have 



T qc (k) 



(20) 



.LAA + _ / dz ^ ( Z ) V (q z , Q, z) i,< (z) 



We can now obtain the result with both intrinsic and ex- 
trinsic contributions in the core electron case from the 
analysis of the intrinsic case by following the derivation 
from Eq. [I?] to Eq. |l9] and replacing - V q (r c ) /uj q by the 
full square parenthesis in Eq. |20j The result is given in 
Eq. |25|. The expression in Eq7^3| for the extrinsic part 
for valence electrons is more involved to simplify. From 
the similar structure of the expressions for the core and 
valence electron cases we can however expect a similar 
qualitative behavior. 



IV. THE USE OF THE ELECTRON ENERGY 
LOSS FUNCTION TO ESTIMATE LOSSES IN 
PHOTOEMISSION 



of the potentials as long as they went to zero at the sur- 
face. A reasonable choice is the Inglcsficld one Q] 

V*{z) = F(z;q„Q)Vf, 

F (z; q z ,Q) = [cos (q z z + <j> q ) - cos ^ q e~ Qz \ 9 (z) , (21) 

which is zero at the surface ( z = 0) and rapidly turns 
into a sine wave. The phase <p q is a certain function of q z 
and Q. This leads us to replace the bulk expression 

JmW(z,z';Q,cj) ~ ^ e ^0" z ') l m W (q z ; Q, uS) 



Qz 



by 



Im W (z, z'; Q, u) 

F (z; q z ,Q)F (z' ; q z ,Q) Im W (q z ; Q, u) . 



For correlated quasi 2D systems we are far from trans- 
lational invariance in the z direction. Suppose we still 
have translational invariance parallel to the planes. We 
write the response between the induced charge density 
and the total (external plus induced) potential as x°> 
P v ' 



— x°V tot , and approximate 



X °(z,z';Q) = X ° b (z~z';Q) 



(22) 

cm — d n ) w (z' — cm — d n ) X2 (Q) • 



It is is a heavy computational task to do a priori cal- 
culations of the fluctuation potentials that are needed 
to evaluate extrinsic and intrinsic losses. Even when 
we limit ourselves (as we have done here) to expres- 
sions quadratic in the potentials, which can be described 
by Im W (r, r'; cj), the task is heavy since it is impor- 
tant to allow for the behavior of W at the surface. To 
avoid these difficulties, awaiting advances in computa- 
tional techniques, we will develop approximations for 
W based on knowledge of the experimental energy loss 
function — Ime^ 1 (q, oj). This is done in three steps. 
First we connect the bulk expression for Im W (r, r'; ui) to 
Ime -1 (q, u>), and then we modify Im W(r, r';ui) to ac- 
count for the surface. Unfortunately not much is known 
about energy loss at finite momenta, and one hence usu- 
ally also have to make the third step of extrapolating the 
experimental Ime^ 1 (q = 0,u>) values to finite q values. 

For a simple metal or a valence semiconductor, an " sp- 
system" (^2|, we take e -1 as diagonal in q space. We then 
have 



iq(r 



>v (q) Ims 1 (q, uj) c?q. 



The fluctuation potentials for a bounded jellium have 
been discussed in Ref. pi] ]. The results for photoemis- 
sion turned out to be fairly insensitive to the precise form 



Here w (z) — \<f>o(z)\ , where <po (z) is the ground state 
wave function for electrons localized in a 2D plane. The 
length of the unit cell in the z direction is c, and d n gives 
the position of the layers inside a unit cell (when there 
are more than one layer). We thus write the response to 
the total potential as one part coming from extended 3D 
excitations plus one part from electrons which are local- 
ized in the z direction and only move in the planes. The 
variable Q comes from the Fourier transform of R — R', 
the coordinate distance in the plane. The susceptibility 
relating the induced charge density to the external poten- 
tial X: P md — xV ext is related to x° by x = X° + X° V X> 
and the inverse dielectric function is e -1 = 1 + u%. Eq. 
p2| leads to a 3D bulk contribution to Im W (z,z'), plus 
a contribution VX2V, where the susceptibility X2 is cal- 
culated as if we only had the 2D planes and then re- 
placing the bare Coulomb potentials in xi with dynam- 
ically screened potentials from the embedding 3D elec- 
trons. We note that the bulk screened potential can be 
anisotropic since x% (Qz> Q) can depend on both q z and 
Q, and not only on q 2 z + Q 2 . If there is only one layer in 
the unit cell there is an explicit solution of the integral 
equation xi = xl + X2 V X 2 > 



X2 (Qz,q z ) = -w (q z ) w (q' z ) x 2 (?*) 
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The variables Q and u> are suppressed, w (q z ) is the 
Fourier transform of w (z) and x 2 (Qz) is a periodic func- 
tion of q z . Clearly w (0) = 1, and often it is reasonable 
to take w (q z ) — 1 for all q z . Due to the periodicity q z 
and q' z must differ by a reciprocal lattice vector. This 
means that if we know the diagonal part of x 2 (<lz,q' z ) 
(from the loss function) we also know the full Xi {izi q' z )- 
With several layers in the cell we have 



X2 (QziQ'z) = \ w (Qz) w ( < l' z 



■E 



Xnn> (<fe) e 



Here Xnn' (q z ) is periodic, but the exponential factor is 
not. The diagonal part of \i (QziQz) contains a sum of 
the matrix elements x nn > (q z ) and thus the loss function 
cannot give us all of them. The diagonal elements of 
Xnn' (<7z) are equal, and if we neglect the non-diagonal el- 
ements (the coupling between the layers in a cell), again 
we have only one unknown quantity, which we can get 
from the loss function. Estimates for Bi2212 show that 
the nondiagonal elements are smaller than the diagonal 
ones but still it is a fairly crude approximation to 
neglect them. In lack of other information than the loss 
function, we however feel forced to do so. 




z (unit:c) 



FIG. 3. The fluctuation potentials V r (z — Zi) and 
V r (z — Z2) where z\ and 22 are the positions of the first and 
second CuO layers in Bi2212. Note that in both cases the 
potentials are zero at the surface, 2 = 0. 

To obtain the 2D contribution to W (z, z') we have to 
sum 

Y / e-^ z e^ z 'v(q z )x 2 (q z , q' z ) v (q' z ) , 

where the variables Q and u> still are suppressed. We 
split the summation over q z into a sum (integral) over 
the Brillouin zone and a sum over reciprocal vectors G, 

W(z,z') - W 3D (z - z') 

,-tt/c 

— /] / V (z — d n , q z ) V (z' - d n ,q z ) xx,i [q z ) dq z , 



where 

V(z,q z )=e- i ^VP(z,q z ), 

VP (z, q z ) = ^ lGz v (q z + G)w (q z + G) . (23) 
G 

The spacial dependence of the fluctuation potential is 
thus that of a Bloch wave, a plane wave times a periodic 
function V p . V p is a well-known function that can be 
obtained explicitly if we take w (q z ) = 1. We define real 
functions 

V r (z, q z ) = Re \ e -^+^,z )y P ^ qz) 

and choose the phaseshift (q z> Zq) to have V r (2, q z ) = 
at the surface z = zq. This procedure is like we would 
have taken F (z; q z ) — sin[q z (z — zq)] which is not unrea- 
sonable. In summary this gives us 

lmW 2D [z,z') 

/7f/c 
V r (z - d n ,q z ) V r (z' - d n ,q z ) Imxi.i (q z )dq z , 

where Imxi.i {qz) can be obtained from the loss func- 
tion. We should have replaced v(q z ) by a screened po- 
tential W (q z ) everywhere, but since we have the same 
replacements in the expression for loss function, this is 
of no consequence. The coupling functions V r ( z) are 
quite different from the functions V q (z) in Eq. |2l| for 
the 3D case. As shown in Fig. 3 they look more like 
a superposition of surface plasmon functions centered at 
the different 2D layers then a sine wave. 

The third step of extrapolating to finite momenta also 
involves considerable uncertainty. It is however impor- 
tant to do some estimates of the effects of dispersion since 
we know that for an electron gas it makes a difference of 
a factor of two for the mean free path. For an electron 
gas the plasmon pole model is fairly good, also when a 
simple quadratic dispersion is used 



Im ■ 



-1 



e(q,u) 



S(UJ - LU q ) , LU q = LU p + q 2 /2. 



Often a one mode description is not good enough. One 
can then use a sum of Drude functions, which further 
include damping of the quasi-boson modes, 



Im 



_L = V C ° 

e(q^) 4-* " 



^r„ (q) 



u n (q) = <J° n + q 2 /2, (24) 



as suggested by Ritchie and Howie Q and others. |5Q] 
The dispersion of oj n (q) is somewhat ad hoc, however for 
large q it has the correct form, following the Bethe ridge. 
Also for r„ (q) some dispersion could be used. The coef- 
ficients c n , and T„ (0) can be determined by fitting 
to the energy loss function — Ime -1 (q = 0, w). 







V. INTERFERENCE BETWEEN INTRINSIC 
AND EXTRINSIC LOSSES 

A. Localized strongly correlated systems 

There is a fundamental difference between photoemis- 
sion from solids and from localized systems like atoms 
or molecules. The interaction between the photoelectron 
and the target always goes to zero at high kinetic ener- 
gies. In solids this decreasing interaction is compensated 
by the longer distance the electron travels, due to its 
longer mean free path, before it leaves the solid (assum- 
ing that the light penetrates much deeper into the solid 
than the mean free path) . Thus while at high energy the 
photoelectron coupling can be neglected for a localized 
system, it never becomes negligible for a solid. 

An extensive literature has appeared on near-threshold 
behavior of satellites in atomic spectra, with a landmark 
paper from 1965 by Carlson and Krause where they 
observed the characteristic shake-up behavior for the first 
time. They showed that the intensity of the shake-up 
satellite in K-shell excitation of the Ne atom rouse slowly 
from zero at threshold to its sudden approximation limit 
at some 200 eV higher energy. The data were plotted in 
units of the " excitation" energy Eq of the satellite final 
state relative to the main line (Eo = 50 eV in the Ne 
case). The onset of satellite intensity was studied par- 
ticularly in the 80's when good synchrotron sources had 
become available. Stohr, Jaeger and Rehr in 1983 |5^] 
made a systematic study of the energy dependence of a 
core level shake-up feature (the K spectrum of N 2 on 
Ni (100)) with a much smaller value of Eq (5 eV). Also 
here the approach to the sudden limit was found to oc- 
cur within a few unit of Eq . Reviews of the atomic work 
on threshold satellites have been given by Becker and 
Shirley. |53|, The picture that evolved was complex 
with rapid variations in satellite intensities near thresh- 
old. Much of the complexity is however related to Ry- 
dberg levels and post collision effects, which have less 
relevance in solids. 

The satellite turn-on was recently studied in detail for a 
numerically exactly soluble model that should be generic 
for a class of strongly correlated solid state systems, jjj 
The model has a charge transfer excitation, and is rele- 
vant for transition metal halides, rare earth compounds, 
chemisorption systems, and high T c compounds. It is 
modelled by three electron levels, one core level and two 
outer levels. When the core hole (3s) is created, the more 
localized outer level (d) is pulled below the less localized 
level (L). The spectrum has a leading peak corresponding 
to a charge transfer between L and d (" shake- down"), and 
a satellite corresponding to no charge transfer. The " con- 
jugate shake-up" mechanism (optical transition between 
two localized levels (s and d) plus shake-up from a local- 
ized level (L) to the continuum) analyzed by Thomas ]55| ] 
is not relevant in our case since an s-d transition is opti- 
cally forbidden. The crucial feature of the model is the 



dynamic Coulomb interaction between the photoelectron 
and the d and L levels, i. e. the interaction that causes 
extrinsic losses. The calculations were done with pa- 
rameters relevant to different copper-dihalide compounds 
(CuBr2, CuCl-i^CuF-i), and showed a rapid approach 
(see Fig. 4) to the sudden limit within 5-10 eV, which 
is in accord with the Stohr et al results. To the extent 
that the localized system described by this model can be 
regarded as decoupled from the rest of the (3D) solid, 
this means that above say 10 eV extrinsic losses and in- 
terference effects are associated only with 3D extended 
excitations. 
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FIG. 4. The ratio r (w) between the satellite and the main 
peak photocurrents for some divalent copper compounds. The 
energy is measured from the threshold for satellite onset, and 
the horizontal lines give the sudden limit. 

With an exactly solvable model some detailed mech- 
anisms for the satellite turn-on could be studied. Since 
the transition takes place at such low energies it turned 
out that the variation of the dipole matrix element gave 
effects comparable to that of the dynamic photoelectron 
coupling. When the dipole matrix effect dominates (in- 
trinsic case), the onset is determined by Eq as long as 
Eo is smaller or comparable to the characteristic energy 
Ed = 1 / R 2 d where Rd gives the length scale of the dipole 
matrix element. For Eq larger than E d it is instead Ed 
which gives the characteristic onset energy. The strength 
of the satellite relative to that of the main line obtained 
with the full model divided by the result when only the 
matrix element was considered, showed curves starting 
at 1 at zero energy which could have sizable humps (ex- 
trinsic effects) at an energy of order Ed- Attempts to 
describe these effects by a semiclassical approximation 
were not successful, and no simple picture in time space 
evolved. Lowest order perturbation theory on the other 
hand succeeded in giving some qualitative understand- 
ing. One conclusion is that Eq may not always be the 
correct scaling parameter for satellite onset. 
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B. Weakly correlated solids 

The transition from the adiabatic to the sudden limit 
in core electron photoemission for weakly correlated sys- 
tems was recently studied. The theory is fully quan- 
tum mechanical and includes the interference between 
intrinsic and extrinsic losses (QM). It is a general theory 
for weakly correlated systems and describes not only 
plasmon losses but losses to density fluctuations in gen- 
eral like those from electron-hole excitations. The actual 
calculations were done for an atom embedded in a semi- 
infinite jellium. For such a simple model explicit and 
detailed results could be obtained. The limit expression 
when the interference term in Eq. pp| was neglected (BS), 
gives essentially the same approximation as discussed by 
Berglund and Spicer |n| and used by Penn. Q Compar- 
ison was also made with results from the semi-classical 
theory (SC). In this theory the photoelectron is taken 
to travel with a steady speed on a straight trajectory, 
and the photoelectron losses are identified with the exci- 
tations caused in the solid by the external perturbation 
from the moving classical photoelectron. In Sect. Ill we 
presented general results, partly based on Ref. H]. 
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FIG. 5. Satellite part of the photoelectron spectrum for Al 
(r s — 2.07) from lowest order perturbation theory showing the 
contribution from bulk plasmons as function of photoelectron 
energy (the zero of the energy scale is at the quasi-particle en- 
ergy), for different photon energies. We show results for three 
different approximations, fully quantum mechanical (QM) 
from Eq. |26| (dotted curves); neglecting interference (BS) 
from Eq. |27| (solid curves); and semiclassical (SC) from Eq. 
pi] (dashed curves). The inset in the ez = 10 curve shows re- 
sults with two different polarization potentials (see Ref. [21]). 

The photocurrent for the core electron case is 



J k (u) 



k|A|c 
1 



dt 
2> 



J{LU-Ek + E c )t 



< ■ " ~ " (25) 



<'Nl> [ - / 9qm (a/) (e 



1) da/ 



We have used Eq. and replaced the Bloch functions 
with plane waves, to find 



9qm M = 
-Vi (z c ) 



E 



- r dze^-^'-^V (z) 



(26) 

5(UJ - UJ q ) . 



Here z c is the position of the atom from which the pho- 
toelectron comes, q labels the fluctuation potentials V q 
(mostly taken from Eq. ||l]), q — (q z , Q), and 

k 2 = k 2 z + 2(\v \ + ir f ), 



= 2 (\V \+w g + ir 2 )-\Qf. 

Further V = (V e ff) average , V eff is the effective crys- 
tal potential, Ff — |ImS (ek + \Vo\)\, and T 2 = 
|ImS(e k + \V \ +u q )\ (c. f. Eq. |). The photoelectron 
is taken to leave the solid at right angle to the surface. 

Satellite spectrum (exponential expression) 

Bechstedt surface and bulk plasmons (ez=50h) 
30.0 , . , . , . , . , . , . , . , 



averaged over core-distance zc 
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FIG. 6. Satellite spectra from the exponential expression 
in Eq. ^ in a plasmon pole approximation neglecting elec- 
tron-hole contributions but including both surface and bulk 
plasmons. The spectra are averaged over core hole position 
and for a high photon energy (ez = 50). 

In the high energy limit the interference terms in Eq. 
p6| drop out and we have |2l| , firj 



9bs (w) = h z c t (e,u) . 



(27) 



where a (to) is the function giving the intrinsic spectrum 
(a (0) is the MND singularity index Q), and r (e, id) is 
the differential inverse mean free path (see e. g. Ref. 
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Eq. g7| gives a result roughly the same as that obtained 
from the transport theory used by Penn. [^9| Finally the 
semi-classical theory with the photoelectron as a perturb- 
ing charge on a straight trajectory gives [glj 



9SC M = 

Vi (z c ) 



E 



+ — / dze-' lu ^ z - z ^l v V q (z) 



(28) 

6(U- UJ q ) , 



where v is the photoelectron velocity. 
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hole pairs and no surface plasmons). Actually the results 
from the exponential expression and lowest order theory 
are identical apart from a constant factor in the energy 
range up to the second satellite peak. The semi classi- 
cal theory works almost perfectly except at lower photon 
energies (see curve with ez = 2) , while the Berglund and 
Spicer transport theory works poorly until the energy is 
very high (of order keV). The quantity ez is the maxi- 
mum energy that the photoelectron can have with a given 
photon energy. 

In Fig. 6 we show the satellite region over an extended 
energy range at a high photon energy (ez = 50), cal- 
culated with the exponential expression in Eq. [2^, and 
averaged over core hole position. We see the strong bulk 
plasmons preceded by the weaker surface plasmons. The 
higher order plasmon peaks become broader and broader. 
The very sharp first and second order peaks will become 
broader when plasmon damping and electron-hole exci- 
tations are accounted for. 
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FIG. 7. Curves of a e ff (to) and P^ (u>) for Al metal with 
a photoelectron energy of 0.85 au (upper panel) and 9.85 au 
(lower panel). The Pk (co) curves are for given k values as 
functions of u, i. e. they are "constant final state spectra" 
(CFS). 
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FIG. 8. The same results as in Fig. 7 for Cu metal. 

In Fig. 5 we show results for the bulk plasmon satel- 
lite in the three approximations. The results are from 
the lowest order term in expanding the exponential in 
Eq. |25| using a plasmon pole approximation (no electron- 



FIG. 9. Curves of the asymmetric quasi-particle peak for 
Al metal (left) and Cu metal (right). The curves have a 
300 meV Lorentzian broadening, and give results for emission 
from an atom 10 au from the surface. Results for a photo- 
electron kinetic energy of 0.85 and 9.85 au, and also when 
only the intrinsic contribution (independent of photoelectron 
energy) is included. 

Al metal is rather special in having a very sharp plas- 
mon peak carrying most of the loss spectral strength (see 
e. g. Fig 4.14 in Ref. §|). To study the effect of the 
shape of the loss function we did some comparisons be- 
tween Al and Cu. We define an "effective" singularity 
function by a e ff (u>) — wgQ m (^) which includes both 
intrinsic and extrinsic amplitudes, and thus depends on 
the photoelectron energy. In Fig. 7 we show a e // (w) for 
Al at two different photoelectron energies, and compare 
results from taking only the intrinsic contribution (intr) 
with the BS and QM cases. Clearly both extrinsic and 
interference effects are strong also at the higher energy. 
In Fig. 7 we also show the photocurrent as function of the 
photon energy obtained with the exponential expression, 
Eq. [ll] at the two photoelectron energies (CFS mode). 
The same results for Cu are shown in Fig. 8. It is clear 
that extrinsic and interference effects are just as strong in 
Cu as in Al. The difference in the loss spectra of the two 
metals however shows in the shape of the satellites, the 
Cu satellite is featureless as could be expected. In Fig. 
9 we show both the intrinsic (intr) and the quantum me- 
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chanical (QM) results for the quasi-particle peaks for Al 
and Cu, broadened by a Lorentzian of FWHM=600 meV, 
and for two photoelectron energies. It is clear that both 
FWHM and asymmetry in practice are unaffected both 
by photoelectron energy and by extrinsic losses. 



C. Strongly correlated quasi two-dimensional solids 

In many high temperature superconductors (HTC) we 
have layered materials with planes of CuO which are ac- 
tive in superconductivity. Much theoretical work con- 
centrates on the electronic correlations in these quasi-2D 
CuO planes, and there are many interesting relations be- 
tween properties found theoretically in purely 2D systems 
and experiment. To interpret photoemission data the full 
3D crystal should however be taken into account. In a re- 
cent study of Bi2212 we find strong effects on the thresh- 
old line shape due to low energy excitations connected 
with the coupling between the planes, effects which do 
not appear in a pure 2D treatment. [Q The effects can 
be seen as external if we consider photoemission from 
one layer and regard the layers outside as an external 
system, or dominantly intrinsic if we look at the system 
as three-dimensional. The effect we discuss is different 
from that considered by Haslinger and Joynt |58[ , who 
consider ohmic losses of the photoelectron occurring af- 
ter the electron has left the solid. 



case, we can however not rest on a comparison with va- 
lence electrons but have to rely on the localization in the 
z direction. 
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FIG. 10. The loss function for Bi2212 as parametrized by 
Norman et al.[59] 

The low energy part of the photoemission spectrum 
comes from exciting electrons in a 2D layer of Bi2212. 
The one electron states are then quite localized in the z 
direction, and the shake-up and loss effects are assumed 
to be similar to those in a core electron spectrum. As 
further motivation for this assumption we can compare 
with results for simple metals. In Sect. Ill we noted that 
the valence electron satellite structure in simple metals is 
very similar to the core electron one. Further the valence 
electron quasi-particle line shape, away from the Fermi 
level, is markedly asymmetric just like a core line. ]46| ] 
For electrons at the Fermi level, which is an important 
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FIG. 11. The contributions to the inverse mean free path 
1/A from the 2D (full drawn) and 3D (dashed) terms in the 
case of Bi2212. 

The calculations for Bi2212 take as input data a three 
term parametrization (see Eq. p4[ ) made by Norman et 
al |59 of the loss spectrum obtained by Niicker et al 
|u0| . The parametrized curve is shown in Fig. 10. The 
first sharp peak describes 2D plasma vibrations, and the 
linear small u> part the acoustic plasmons coming from 
coupling between the 2D layers. In Fig. 11 we show the 
contributions to inverse mean free path coming from the 
2D part (the first term), and the 3D part (the second 
and third terms). It is clear that at about 20 eV, where 
photoemission measurements usually are done, the 3D 
contribution dominates. As discussed in the introduc- 
tion, we take the statevectors as direct products, and 
write for the initial and final states, 



\N B ) \N 2D ) , \N* B , Sl ) \N 2D - 1, s 2 ) |k> . 



(29) 



Here \N2d) is the state vector for the electrons in one 
particular layer at a distance zq (zq > 0) from the sur- 
face (the one from which the photoelectron comes), and 
\Nb) the state vector for the remaining (bulk) electrons 
which move in 3D. \N 2 d — 1, s 2 ) is an excited state s 2 of 
the particular layer, and \N B ,Si) an excited state si of 
the bulk electrons. The star indicates that these electrons 
move in the presence of a localized hole at r = (0, zo), 
and thus is an cigenfunction of a different Hamiltonian 
than that for \Nb)- Finally |k) is the photoelectron state. 

Assuming that the 2D part is like a core level, we could 
derive an expression for the photocurrent as a convolu- 
tion between the 2D current J^ D (zq,ui) and an effective 
broadening function (zo, w), 

Jk {zo, uJ phot ) = J Jl D (z , lu') P k (z , Lo phot - w') dJ . 

(30) 
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A delta function peak 5 (lu — eo — £k) in J%P (z , ui) will 
hence give a contribution Pk (zQ,Lu p hot — £o — £k) to the 
photocurrent. Summing up to an exponential expression 
we have 



(31) 



x exp 



a 2D (k, 2o;w') ■ 



-duj' 



"3D (k, z ;w') ■ 



-dui' 



The 3D contribution is somewhat arbitrarily cut at a 
threshold value of to t h =0.1. We have replaced e 



iu t 



1 



by e tuJ * in the 3D term. The reason is that in the high 
energy limit we have 



"3D (k, z ; oj) 



dtu = 



zq/X, 



(32) 



provided we neglect the contribution to the mean free 
path A coming from 2D excitations. This, as we just re- 
marked, is reasonable at photoelectron energies of 20 eV 
and higher. Since the 2D and 3D contributions add in 
an exponent we can write Pk (zq,lu) as a convolution, 



PkC^oX 



p2D 



For P^ D we make a Taylor expansion, and keep only the 
low order result, P^ D (z , oj) = 5 (uj)+a 3 £) (z , lu) /lu. We 
have then omitted the multiple quasi-boson excitations 
starting at lu = lioth- Since P£ D is normalized to unity, 
and consists of a peak that is sharp compared to «3D, wc 
can write 



-Pk (z ,u) 



-z /A- 



p2D 



a 3D (k, z ;uj) 



(33) 



The extrinsic contribution to a 2 o (k, Zq)uj) can be ne- 
glected, and the intrinsic contribution has a zq but no 
k dependence. There is a competition between a de- 
clining function exp(— zq/X) and an increasing function 
(%2d{zo,uj) making the contribution from the second 
layer dominate. We approximate a 2 o (zq,lu) by a rect- 
angular function, aQ0 (luq — lu) independent of zq, and 
broaden with a Lorentzian of width T (FWHM = 2T). 
As parameters we use ao — 0.255 and luq = 0.08. For 
lu < ujq we then have the Doniach-Sunjic expression [ pi[ , 



P k 2C H =C(oo) 



C(a ) 



cos [7rao/2 — (1 — a ) arctan (w/T)] 



(l + W 



1 

-ya 



(1-q )/2 



(a - ljlWoT 1 -"" sin[7ra ] 



(34) 



where 7 = 0.577 is the Euler constant. For lu > luq 
P^ (a;) only has a weak tail with less than 10% of the 
norm (for ao < 0.4). The coefficient C (ao) was derived 
in Rcf. ii. 
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FIG. 12. The effective loss function P^ D (lj) convoluted 
with Lorentzians of different widths V (5, 10, and 20 meV). 
Also the Lorentzians are shown. The photoelectron energy is 
0.85 au. 



P£ D (lu) is the function that broadens a ^-function 
peak in j£° (lu). If J^P (lu) has a Doniach-Sunjic sin- 
gular shape the broadening with P£ D (lu) still gives Eq. 
BJ but with an ao that is the sum of the alphas in j£ D 
and in P£ D (lu). This is so because the time transform of 
a power law singularity t a ° , and a convolu- 

tion in frequency space is a product in time space. 

In Fig. 12 we show the sum for the first four layers of 
the 2D contributions exp (— z /X - a 3 n ^ r ) P£ D (lu) broad- 
ened with different Lorentzians. The 3D terms are not 
included except for the (all important) mean free path 
factor. In the first three panels with a limited energy 
region (up to 500 meV) we have used the rectangular ap- 
proximation for the a 2 D contributions. In the last panel 
with a larger energy range the full evaluation was done 
since it is superior to the rectangular model for energies 
above 0.5 eV. Also the Lorentzians are shown to ease the 
estimate of the size of the asymmetries. It is clear that 
we have a sizeable line asymmetry, and also a long tail 
extending over several eV. 
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In the superconducting state the loss function should 
have a gap. We mimic this gap by using a rectangular 
alpha function 

OL2D (y>) = a o 0(cj - lo sc )6 (ui - w) , (35) 

still using a = 0.255 and lo = 0.08 au = 2.2 eV. For 
the gap u sc we take lu sc = 70 meV. In Fig. 13 we show 
the corresponding P£ D (oj) broaden with a Lorentzian 
of width r = 15 meV. Our choice of parameters is 
only made to illustrate the qualitative behavior to be 
expected. The curve clearly shows the peak-dip-hump 
lineshape found experimentally (for a recent reference see 
eg ©) 
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FIG. 13. The effective loss function P% D (w) for a gapped 
spectrum using the parametrization in Eq. The 
Lorentzian broadening is V = 15 meV. 

Recently it has been possible to obtain very accurate 
tunneling data from Bi2212, and it is of interest to com- 
pare these data with the PES satellites (Ref. |6J|), since 
the tunneling data also show peak-dip- hump structures. 
Ipij PES and tunneling are basically different spectro- 
scopies. There can however be qualitative similarities 
since in both cases the electrons couple to 3D quasi-boson 
excitations like phonons, electron-hole pairs, plasmons, 
magnons etc. In our treatment of PES we take the prop- 
erties of a particular 2D layer as given and study the 
effect to low order of the sudden appearance of a hole in 
the 2D system on the quasi-bosons (intrinsic excitations) 
as well as of the coupling of the photoelectron leaving 
this layer to the quasi-bosons (extrinsic excitations), and 
their interference. We found that the intrinsic contribu- 
tions dominate for small excitation energies. 

Tunneling is traditionally described by a spectral func- 
tion which involves matrix elements of the electron anni- 
hilation operator between the initial state and the excited 
states. [^3| The excited states consist of a 2D layer 
state with a hole, and some state of the quasi-bosons 
in the presence of a localized hole. In lowest order per- 
turbation theory the probability for a final state with 
excited quasi-bosons is given by the first term in Eq. |2(]. 
This means that the intrinsic contribution to PES and 
the tunneling currents are the same except for the mean 



free path effect in PES shown in Eq. ^3[ and the sum- 
mation over momenta in tunneling giving the density of 
states, DOS. In PES on the other hand, we have momen- 
tum conservation from the dipole matrix element. As 
mentioned above, we modify this analysis valid for the 
normal state, by simply assuming that the loss function 
should have a gap in the superconducting state. 

In Bz2212 we have a van Hove singularity (VHS) at the 
Fermi level, which makes the difference between DOS and 
momentum conservation of less importance (there might 
actually even be two VHS if the two CuO planes at 3 A 
apart produce a significant splitting). More important 
is that in PES the electrons come from a thin surface 
region (of the order of the mean free path) while in tun- 
neling they may come from an extended region which 
can be hundreds of A, and that the coupling functions 
V(z) have a slow approach to their bulk limit (Fig. 3). 
Additionally there are two energy gaps (superconduct- 
ing gap and pseudogap), which further complicates the 
picture. It is clear that PES and tunneling structures 
cannot be quantitatively the same, but since the same 
quasi-bosons are involved, there may well be qualitative 
similarities even though the coupling strengths can be 
different. We remind that we have taken the spectral 
function for the 2D system (the function often calculated 
by theoreticians using say a t-J model) as a sharp peak, 
and any peak structure in Jid (w) has to be convoluted 
with the effective broadening function P^d ■ In our anal- 
ysis we have only treated the plasmons for the simple 
reason that the experimental loss data at hand did not 
have resolution enough to show phonons and other low 
energy excitations. 



VI. CONCLUDING REMARKS 

When we go beyond the common procedure of com- 
paring the spectral function with photoemission, there 
are a number of definitions and distinctions that have to 
be introduced. In the one-step approximation only the 
quasi-particle peak and the dipole matrix element are in- 
cluded but no intrinsic or extrinsic losses. In the intrinsic 
approximation the full spectral function including its tail 
from intrinsic losses, as well as the dipole matrix ele- 
ment are included, but no extrinsic losses. In the sudden 
approximation (or Berglund-Spicer model) the spectral 
function is taken as a source term in a transport prob- 
lem, or roughly equivalently we make a convolution be- 
tween spectral function and the extrinsic loss function. 
The dipole matrix element is taken as independent of 
loss energy. The adiabatic approximation somehow de- 
scribes the behavior at threshold, but there is no simple 
approximation scheme to do this. When the losses come 
from extended excitations the satellite in photoemission 
has contributions both from intrinsic and extrinsic am- 
plitudes. A semi- classical approximation puts the pho- 
toelectron on a trajectory, and the losses are calculated 
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with the photoelectron as an external perturbation. It 
should be noted that there is considerable ambiguity in 
the literature concerning these definitions. 

In this paper we discuss photoemission in general, we 
explore the limits of the sudden approximation, and we 
analyze the role of extrinsic losses. It is not a compre- 
hensive treatment, that would require much more space, 
and we concentrate on some problems of broad interest. 
Thus in core electron photoemission we e.g. do not take 
account of the finite core electron lifetime. This lifetime 
can cause spectacular effects if phonon energies are of the 
same order as the lifetime width. It can also cause 
large effects if there e.g. are Koster-Kronig decay chan- 
nels, and two nearby discrete levels (say 2pi/ 2 and 2p 3 / 2 )- 
|37j There are many resonace effects in d and f electron 
systems, some involving atoms on different sites like mul- 
tiatom resonance PES. [H We have omitted this large 
topic. For a general review see Ref. ]69| ] 

We study in some detail how the three-dimensionality 
of a quasi-2D solid, Bi2212, influences its photoemission. 
We have earlier found that for a strongly correlated local- 
ized system the sudden approximation is reached rather 
quickly, at about 10 eV ffl. For a weakly correlated 
system on the other hand, like an sp-metal or semicon- 
ductor, the sudden limit in core electron photoemission 
is approached very slowly, on the keV scale. M] The 
slow approach is connected with strong destructive in- 
terference between the intrinsic and extrinsic channels 
for plasmon production. The cancellation is particularly 
strong for small momentum plasmons where the long- 
wave plasmons are excited by the average potential from 
the core hole and photoelectron, which is zero [^0| . When 
we compare the satellites of Al and Cu we find that the 
latter have about the same strength but are structureless. 
The asymmetric lineshape in core electron photoemission 
from metals is, on the other hand, hardly affected by the 
external loss processes plj| . 

We are interested in energies where the sudden limit 
is reached for the strongly correlated layer from which 
the photoelectron comes, and derive an expression for 
the photocurrent as a convolution of the sudden approx- 
imation for the current from the layer with an effective 
loss function, Pk (w) (Eq. ||(]). We assume, as far as 
the loss properties are concerned, that the photoelectron 
comes from a localized position. In our specific example, 
Bi2212, the c value is 15.4 A (neglecting crystallographic 
shear), and almost all contributions come from the first 
unit cell. The two first CuO layers are at 0.39 c and 
0.61 c from the surface (which is between two BiO lay- 
crs). @ _ _ 

To obtain Pk (w) we use a previously developed method 
based on a quasi-boson model, where the electron-boson 
coupling is given by fluctuation potentials related to the 
dielectric response function |^] . We find that the energy 
loss function, which we take from experimental data, can 
be related to the screened potential which we need to 
calculate the (intrinsic and extrinsic) losses in photoemis- 
sion. The fluctuation potentials related to the electrons 



in the layers are universal functions, which are easily cal- 
culated. They have some resemblance to a surface plas- 
mon potential, but penetrate the whole solid and have the 
Bloch wave symmetry. We use the real part (or equiva- 
lently the imaginary part) of a phase shifted bulk poten- 
tial to get a potential which is zero at the surface, and 
mimics the potential we have in a finite solid. 

We have extrapolated results from electron energy loss 
experiments for q=0 to finite vales of q. This inclusion of 
dispersion makes the mean free path considerably longer 
than obtained by Norman et al j5j|, about 12 A rather 
than 3 A, at say 20 eV (Fig. 2). Measurements by the 
ITR-2PP technique Q give a lifetime of r = 10 fs at 
an energy e = 3 eV above the Fermi surface. The mean 
free path is A = vt. Converting energy to velocity by 
mv 2 /2 = e gives a mean free path r = 103 A as com- 
pared to our result of about 17 A at that energy. This is 
an indication that our values rather are on the low side. 
It is however hard to know what is the correct conversion 
between energy and velocity at such low energies, which 
makes a comparison very uncertain. 

From Fig. 11 we see that the 2D losses occur only for 
small energies, at 5 eV the bulk losses take over. The 2D 
losses go to zero quite slowly, just like the bulk losses, 
but on another energy scale. If we only had 2D losses, 
the minimum mean free path would be long, about 20 A. 
The general behavior of the 3D mean free path follows 
a well known pattern. The mean free path has a mini- 
mum of about 5 A at an energy of 3-4 times the energy 
where the loss function has its center of gravity. We 
have used the Born approximation to evaluate the mean 
free paths. This may seem a very crude approximation 
at low energies. However the Born scattering expression 
with a basis of Bloch waves and Bloch energies rather 
than plane waves and free electron energies agrees with 
the GW approximation, which is commonly used also at 
low energies. Further it was shown by Campillo et al |73j | 
that plane waves and free electron energies was not that 
bad, as long as the energies in the dielectric function are 
well approximated. 

Our main concern is the behavior of the effective 
broadening function at small energies where it is dom- 
inated by the 2D losses. The 3D contributions set in at 
somewhat higher energies, and give a rather structure- 
less contribution. What we here for convenience call 2D 
losses is of course actually also a 3D effect since it comes 
from excitations of a coupled set of 2D layers. To allow 
a qualitative discussion we represent the a 2 D functions 
by a rectangular distribution. The rectangular distribu- 
tion allows an analytic solution valid out to the cut-off 
ujo (u>o ~ 0.08 au ~ 2.2 eV). Pk(u>) has only a fairly 
small tail beyond u>o. In Fig. 12 we plot the total Pk (ui) 
function (sum over the four first layers, properly mean 
free path weighted), calculated with the rectangular ap- 
proximation and broadened with Lorentzians of different 
widths. We note the marked asymmetry. If the J 2 d 
function has a power law singularity with singularity in- 
dex aL, the photoemission will have a singularity index 
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(a + a L ). 

In a paper by Liu, Anderson and Allen from 1991 [f74| , 
they discussed the lineshapes of BiiSriBaCuiOz along 
the T — X direction obtained by Olsen et al for 22 eV 
photons. They concluded that neither the Fermi liquid 
nor the marginal Fermi liquid theories could fit the slow 
fall-off of the spectrum at higher energies. Our results 
offer a possibility that the slow falloff may be due to in- 
trinsic creation of acoustic plasmons in a coupled set of 
CuO layers, an effect not present if only one CuO layer 
is considered. This broadening is mostly intrinsic, i e if 
we treat a 3D system we have an almost intrinsic effect. 
However most theoretical discussions concern an isolated 
2D system, compared to which we find an appreciable 
extra broadening from the coupling between the layers. 

The PES spectra change strongly when we go to the 
superconducting state. The main peak sharpens and a 
peak-dip-hump structure develops. This effect has been 
interpreted as a coupling of the 2D state to the (tt, it) 
collective mode. ]6^] Here we find that this effect also 
can arise from the gapping of the loss function caused by 
the lack of low energy excitations in a superconductor as 
shown in Fig. 13. Without a more accurate model we 
find it difficult to decide which is the correct explanation, 
possibly it could be a combination of both mechanisms. 
Since the gapping of the loss function is related to the su- 
perconducting gap, also with our mechanism the hump 
will scale with the gap. It is clear that the experimental 
peak-dip-hump structure rides on a background which is 
not predicted by our expressions, nor by anyone else's. 
Our theory is however rather schematic with its strict 
separation of a 2D and a 3D part, while in reality the 
bands are hybridized. If we extend our approach to a 
more detailed treatment of the underlying bandstructure, 
the background could well be strongly changed. Such an 
extension represents a very large numerical task, but with 
the present pilot treatment we can at least start think- 
ing seriously about the difficult background problems in 
photoemission. 

In a recent paper Haslinger and Joynt |5£| discussed 
a broadening mechanism due to the interaction between 
the photoelectron when outside the solid and the elec- 
trons in the solid. This is a different mechanism than 
in this paper, which adds additional broadening. This 
mechanism has also been discussed by Schulte et. al. |76| ] 

It should be stressed that we cannot claim any high 
quantitative accuracy. We have put in dispersion in the 
loss function using a crude approximation. Since how- 
ever dispersion is very important we think our predictions 
are substantially better than if dispersion had been ne- 
glected. We have only considered normal emission where 
the electrons come from the T point, while the interest- 
ing experiments concern electrons from the Fermi surface. 
However there is no reason that the effective loss func- 
tion should change qualitatively when we go away from 
normal emission. The behavior of the loss function when 
lo — > has been disputed. Most authors seem to believe 
the approach is linear, but there are also claims that it 



should be quadratic. 77 If it were quadratic, the cor- 
responding a-function would start linearly rather than 
with a constant. However a (to) would have to rise very 
fast to reproduce the behavior of the loss function for 
the (quite small) energies where it is known to be ap- 
proximately linear. Thus the pure power law behavior of 
Pfc (to) would be lost, but Lorentzian broadened curves 
would probably not differ much. Our fluctuation poten- 
tials are obtained by phase-shifting bulk potentials to 
make them zero at the surface, and define them as zero 
outside the solid. This procedure turned out to be fairly 
good in the metallic case, where we could check with 
more accurately calculated fluctuation potentials. Again 
this approximation is crude, but we believe it to be fun- 
damentally better than if we had used a step function on 
the bulk potential. Since the phase of the bulk potential 
is arbitrary, such a procedure would anyhow have been 
arbitrary. To calculate more accurate potentials is a very 
large numerical undertaking. 

One may also question the use of a bulk expression to 
estimate of the mean free path at the fairly low energies 
that we are concerned with, after all we found strong ef- 
fects when modifying the fluctuation potentials for sur- 
face effects. It is not easy to make a strong statement 
here, and we can only refer to "the state of the art" , that 
bulk mean free paths are successfully used in LEED and 
also in low energy life time calculations which are com- 
pared with time-resolved two-photon PES (TR-2PPE) 
experiments. ]7S[ ] 
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VII. APPENDIX A 

The quasi-boson parameters are related to the dynam- 
ically screened potential W(r, r';uj). Using its spectral 
resolution we can write W in terms of the fluctuation 
potentials V q , 



W{r,r';w) 



e 1 {y,y";uj) 



dv" 



(36) 



9 



2uj q V q (r) V q (r') 



2 2 
0J Z — Li)* 



where, 



y? ( r H / nr^, p*(r) = (N,q\p op (r)\N), (37) 



Pop (r) = / ^ {x)tl> (x) de, u q = E(N, g) - E(N, 0) - 18, 
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and x stands for three space coordinates r, and £ for a 
spin coordinate. Since u q and V q are not known exactly 
we use the RPA where 



V q (v) = / W(T,T f ;u> 9 )4 q (/)dr f , 



(38) 



with uj q = e k -ei (e k > n > ej) and (f) q (r) = cj> k (r) 0; (r) , 
i e the state q corresponds to a particle-hole excitation. 
Undamped plasmons represent a singular case, which re- 
quires a special treatment. |^H , |22| In our treatment 
we will take the state "q" as a quasi-boson state. In our 
applications to photoemission the fluctuation potentials 
V q always appear in a quadratic form together with to q 
such that they can be eliminated and replaced with the 
imaginary part of the screened potential Im W 

ImW (r,r';uj) = -tt E V q (r) V q (r') 5 (lu - u q ) , u > 0. 
g 

(39) 

Im W in turn can be estimated from experimental elec- 
tron energy loss data. 

VIII. APPENDIX B 

To find the ground state for the Hamiltonian in Eq. [l^ 
we replace c\c v by Sw — c v c\ to have 

H syst = E e A c i + E u q a\a q (40) 

l q 

- E v?v<h4 K + 4] + E v u K + <\ ■ 

qlV ql 

The last term can be combined with the quadratic boson 
term, and the linear boson operator eliminated by a shift, 

Vr 



= a - -12- V q - V V q 



The ground state for H syst is thus a filled Fermi sea with 
no quasi-bosons a q , i. e. a simple Slatcrdcterminant 

|iV /p ), and its energy is E a = £, % + E ? W) 2 Av 
This is also the ground state for H since it does not 
contain any photoelectron and thus both h and V give 
zero. The transformation matrix that diagonalizes the 
c|c ; term deviates from 6w by terms of order (V 9 ) 2 , and 
we will neglect the difference between V q v and V q v and 
between ei and ej. When the occupied state I is a core 
state the transformation matrix is unity and V, q , = V, q , . 

The state q \N ip } appearing in Eq. [l^ is 
not an eigenstate of H syst . We note that states 

{d\N IP )}, {a\ Cl \N IP )}, {Z\a\, ci \N IP )) form a 

complete set to span \N — 1, s) . We have to lowest or- 
der in the coupling functions V q 



10 = 10° + ^ 



v q 

1 rv 



^ >'> q - ei> + £( 



W) 



qV 



where 



\l? = c l \N IP ), \ql)° =a\ Cl \N IP ) 



We express the "initial state" \T) in eigenstates of H syst , 



io° = io-E 



v 



ql 



1 u q - e v + e x 



W) 



We need the matrix elements T s (c.f. Eq. fLa), 



T„ = ( N-1,8 



1 + V 



E - H syst -h- QVQ 



N 



ip 



We are interested in the two cases when " a" corresponds 
to \l) and to \ql), 



\N-l,s) = \l) 
E(N-l,s)=E(N)- £l 



to obtain 

H syst = E £lC i c i + E W 9^9 



(41) 



qlV 



~ ~t „v q 

Ulq 



E 



The last term in the square bracket is combined with the 
first term, which then is diagonalized. This gives 



Hsyst = E e A C l + E W 9°9 a 9 

I q 



(V T 



|JV-l,«) = |gi) 
£? (JV - 1, s) = B (AT) - e, + w, 



We have, neglecting O (V 2 ) terms, 

1 



T q l = ( ?Z 



1 + V 



£ - Hsyst -h- QVQ 



I') =6 W 



1 + V 



1 



E-Hs 



syst - h - QVQ 
-V q 1 



uj q - si+ £ V "" Lu q +£u- h - S (w g + e k ) 
Combining these results with Eqs|l| and [l^ we have, 
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J u (w) = ^ |r ; (k)| 2 S (uj - e k + E/ ) (42) 
+ ^ |T g ; (k)| 2 <5 (w - uj q - e k + ej) 

ql 

with 

7i(k) = (k|A|l) 

(k) = £ V 7 (43) 

+ Vq T^TT rA A . 

\ uj q + e k - h - E (u; g + e k ) / 

Here the first term is the intrinsic and the second the 
extrinsic amplitude. For photoemission from a core level 
only one state is involved, and we can replace all I by 
one single index, say c. Our expression then reduces to 
our earlier result in Refs. [^(J and Since r Q ; (k) is 

quadratic in the fluctuation potentials V q Eq. O can 
also be written in terms of Im W . 
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